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ABSTRACT 

Aims. This is the first paper in a series that deals with solar-physics applications of the equation-of-state formalism based on the 
formulation of the so-called "Feynman-Kac (FK) representation". Here, the FK equation of state is presented and adapted for solar 
applications. Its domain of validity is assessed. The practical application to the Sun will be dealt with in Paper II. Paper III will extend 
the current FK formalism to a higher order. 

Methods. A recent rigorous quantum-statistical formalism for Coulomb systems is used to compute the thermodynamical quantities 
for solar modeling, taking into account the necessary requirements on smoothness and accuracy. The FK formalism being a virial 
expansion, it suffers from the well-known deficiency that it is limited to nearly full ionization. This point is elaborated in detail, and 
the quantitative criterion for the domain of validity of the FK equation of state is established. 

Results. Use of the FK equation of state is limited to physical conditions for which more than 90% of helium is ionized. This incudes 
the inner region of the Sun out to about .98 of the solar radius. Despite this limitation, in the parts of the Sun where it is applicable, the 
FK equation of state has the power to be more accurate than the equations of state currently used in solar modeling. The FK approach 
is especially suited to study physical effects such as Coulomb screening, bound states, the onset of recombination of fully ionized 
species, as well as diffraction and exchange effects. 

Conclusions. Despite technical difficulties in its application, there are unique features in the FK approach that promise to turn it into 
the most exact of the available formalisms, provided FK is restricted to the deeper layers of the Sun where more than 90% of helium 
is ionized. The localizing power of helioseismology allows a test of the FK equation of state. Such a test will be beneficial both for 
better solar models and for tighter solar constraints of the equation of state. 

Key words. Equation of State; Sun: interior; Sun: helioseismology 
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, 1. Introduction consistent, valid over a large range of temperatures and densi- 
ties, and incorporate the most important astrophysically relevant 
This paper is the first of a series of three devoted to the study of chemical elements. Let us briefly review the two major equation- 
the equation of state based on the formalism of the path integral of . state efforts of the last 20 years. Both were made in the con- 
in the framework of the Feynman-Kac (hereinafter FK) represen- text of the two most recent opacity recalculations, 
tation. This formalism leads to a virial expansion of the thermo- Gne of me m is the international Opacity Project (OP; see 
dynamic functions in powersofthe total density of^oulomb me books by ge^nl 1995; Berringto^ |S); it contains 
plasma ([Alastuey & Perezj |L992t | Alastuey et all 1199J [19951 the so-called Miha las-Humme r-Dappen equation of state 
l Alastuey&Perezt |1996|). The path-integral approach is the lat- g7j^ mer & Mihalas, 1988; Mihalas, Dappen & Hummerl 



est in a series of attempts to improve the equation of state for 1988; Dap pe n etali fpMH iNavfonov et all [19991: 



the solar interiors. However, so far there has only been a very | Trampe dach et all I2006L hereinafter MHD), dealing with 

limited amount of astrophysical applications of that formalism. heuristic concep ts about the modification of atoms and ions in 

One ( iPerez&Chabneii |2004|) was an application to the screen- a plasma SpecincallVi in MHD the destruction of atomic states 

ing enhancement ofjiuclear fusion. Another one was a precursor is modeledi on the one hand by neutral species according to 

d Perez & Dapper* [1995|) of the present series of articles, which excluded volume> and on the other hand by Stark disso i ution due 

addressed the possibility of using the FK formalism in solar and tQ the surrounding microfie ld. From the included 15 astrophysi- 

stellar modeling. Here, we embark on a more systematic apphca- caU releyant elements (H He c Nj Q Ne? Na? Mg> A1> Si> s> 

tion of the FK representation. It is well known that solar and stel- M ^ Ca> Fe)> users can select any subset Although MHD in its 

lar modeling requires thermodynamic quantities that are smooth, current yersion does not contain free parameterSi several of the 

physical quantities are, in principle, open to revision. Possible 

Send offprint requests to: Asher Perez candidates for adjustment are (i) the atomic and ionic radii 
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(which can be chosen to match observations), (ii) the strength 
of the assumed microfield distribution, and (iii) the specific 
form of the interaction potential between species. Since MHD is 
basically a heuristic equation of state, such future improvements 
will be legitimate and they can be de veloped under the guidanc e 
of observations. For more details see Trampedac h et al.\ d2006). 

The other recent opacity effort is the OPAL project pursued 
at Livermore. In contrast to OP, its underlying equation of state 
has relied on the so-called physical picture, which is built on the 
notion of fundamental particles only, that is, electrons and nu- 
clei. It provides a systematic method to include nonideal effects. 
Its s o-called ACTEX ("activity e x pansion") equation of state 
(see | Rogersl Il986t iRogers et all 119961; iRogers & Navfonovl 



2002, and references th erein) became part of the OPAL opac- 
ity project a t Livermore dlglesias et al. esias & Ro gers. 
1991. 119931 fl995lll996l) . For this reason, the ACTEX equation 
of state is hereinafter referred to as the OPAL equation of state. 
As mentioned, OPAL is a physical-picture formalism, starting 
out from the grand canonical ensemble of a system of the ba- 
sic constituents (electrons and nuclei) interacting through the 
Coulomb potential. In the currently released version, the OPAL 
equation of state includes the following astrophysically relevant 
elements: H, He, C, N, O, Ne (but the OPAL opacity calcula- 
tion includes other elements, such as Fe, as well). In OPAL, 
any effects of the plasma environment on the internal states are 
obtained directly from the statistical-mechanical analysis, rather 
than by assertion as in the chemical picture. Therefore, in con- 
trast to the intuitive MHD, in OPAL there are no adjustable pa- 
rameters. On the one hand, this is its strength, on the other hand, 
it means that one cannot easily upgrade OPAL even if observa- 
tions were to suggest so. Thus in the current OPAL equation of 
state tables, any residual inaccuracy is either structural, a result 
of (i) the finite truncation in the underlying activity expansion, 
and (ii) the utilization of effective atomic and ionic potentials, or 
it could be due to the limited number of chemical elements in 
the current tables. 

The equation of state for solar and stellar structure and 
seismology has to be formally precise and consistent, even 
before the question of the accuracy of the physical descrip- 
tion is asked. It has to satisfy four conditions: i) a large 
domain of applicability (in p, T), ii) a high precision of 
its numerical realization, iii) consistency between the thermo- 
dynamic quantities, and iv) the possibility to take into ac- 
count relatively complex mixtures with at least several of 
the more abundant chemical elements. These requirement are 
generally easy to satisfy with relatively simple formalisms, 
such as the popular, widely-used CE FF equation of state 
dChristensen-Dalsgaard & Dappenl[l992l) . which is an offspring 
of the wid ely used Eggleton-Faulkner-Flannery ( EFF) equation 
of state dEggleton. Faulkner & Flannerv L 1 19731) , or the other 
offspr ing of EFF, the SIREFF equation of state (IRogers et ali 
119961) . 

However, the demands from helioseismology are such that 
formal considerations are not sufficient. The equation of state 
has to be accurate, not just smooth. Even state-of-the-art 
equations of state such as OPA L and MHD are not com- 
pletely sufficient dPappenl 120061) . and improvements are still 
necessary. Improvements can be made either formally (phe- 
nomenologically) or rigorously. A form al option is to merge 
the good feature s of OPAL and MH D dStarostin et ali 120031: 
iLiang & Dapperl 120031; iLiangl l2004t iDappenl 120091) to create 
an equation of state that matches reality as closely as possible. 
However, as a phenomenological effort, such a hybrid formalism 
would not be a scientific improvement but merely a better com- 



putational tool for solar and stellar modeling. As far as the de- 
velopment of a real theory is concerned, genuine improvements 
in the physical model must be made. The rigorous equation of 
state in the FK representation is such an improvement. In prin- 
ciple, it is a better theory than either MHD or OPAL. A case in 
point is the treatment of Coulomb screening. Here, the electrons 
are taken into account with the exact quantum-mechanical ex- 
change and diffraction terms. These are corrections that deal, on 
the one hand, with the "Pauli blocking" of electrons (exchange 
terms) and, on the other hand, their finite thermal de Broglie 
wavelength (diffraction terms). For these terms, both MHD and 
OPAL use approximations. However, this better physics of FK 
comes with a price: the FK equation of state has a limited do- 
main of validity unlike the more global MHD and OPAL. It only 
works for nearly fully ionized plasmas. 



As a theory, FK is quite similar to OPAL, since it is real- 
ized in the physical picture and the grand canonical ensemble. 
However, its net result is that of a virial (density) expansion (un- 
like OPAL, which is based on activity). The virial coefficients of 
the FK expansion are evaluated by the path-integral formalism, 
which in its essence uses the equivalence between a point-charge 
quantum system and a classical one made of extended filaments. 
In this spirit, the calculation of all non-ideal contributions to the 
equation of state becomes possible, systematically, exactly, and 
analytically. The fact that FK has the form of a virial expansion is 
the source of its limitation to near full ionization. In Section 2 we 
illustrate the well-known fact that in the case of reacting gases, 
virial expansions are inadequate to deal with recombination (un- 
like the activity expansions of OPAL). We discuss this deficiency 
quite in detail and develop a quantitative criterion for the valid- 
ity of the FK equation of state. The adopted criterion is a limita- 
tion to physical conditions such that more than 90% of helium is 
ionized. The subsequent sections discuss the physical condition 
quantitatively, in terms of a solar model. 



Despite the aforementioned limitation to nearly fully-ionized 
matter, the FK approach can be applied to solar modeling. The 
reason for this is twofold. First, a large part of the solar interior 
(out to about .98 of the solar radius) lies inside its domain of 
validity, and second, helioseismology allows a localized analy- 
sis of thermodynamic properties despite the fact that there is no 
complete coverage of the Sun by the FK equation of state. 



Here the FK equation of state is presented and its domain of 
validity discussed. Paper II will then be dedicated to solar appli- 
cations. They will be based on "mixed" thermodynamic tables, 
with FK in the highly-ionized part, and a conventional equation 
of state (OPAL or MHD) outside the domain of validity of FK. 
We note in passing that once these tables exist they can be used 
also for stellar modeling, but it is not clear if stellar application 
would ever require such an accuracy of the equation of state, in 
contrast to solar appli cation for whic h this need has been well 
established (see, e.g., Dappen, 2006). In any case, stellar ap- 
plications are beyond the scope of this current series of papers. 
Finally, paper III will be a natural extension of the current level 
of the FK formalism. While the current level of FK is based on 
an expansion of pressure in terms of the density with up to order 
p 5 ! 2 , the next (rather elaborate) step is the extension to calcu- 
lation of the thermodynamic functions up to order p 3 . Then we 
will be able to take into account three-body effects, allowing us 
to treat helium exactly, beyond the hydrogenic approximation. 
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2. Illustration of the limitation of virial expansions 
for reacting systems 

The principle can be understood by considering the exact equa- 
tion of the ideal (non-relativistic) Fermi gas. One starts out from 
an implicit equation of stat e, obtained f rom the grand-canonical 
partition function (see, e.g. Huangl [l963l) 
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Here, p, T, N, V, k, A are pressure, temperature, particle number, 
volume, Boltzmann constant, and the thermal de Broglie length 
A = h/ y/2nmkT, respectively. The fugacity z = exp(p/(kT)) (p. 
being the chemical potential) is defined implicitly by Eq. (1), by 
virtue of the Fermi integrals 
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Expanding the Fermi integrals as a series in the high-temperature 
limit, and eliminating the fugacity 7, one obtains the virial (den- 
sity) expansion 
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(4) 



Alternatively, one can continue to work with fugacity expansion 
1 
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With the introduction of the activity 7 = z/A 3 we obtain an even 
more convenient expansion 
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These expressions above are exact for the partially degenerate 
ideal Fermi gas, they retain their formal structure for any real 
gas (this follows from their origin in the grand-canonical par- 
tition function). Therefore all the physics lies in the value of 
the coefficients of the expansions. For the activity, one formally 
writes 



P_ 

kT 
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7 - b 2 z + 



= z- 2b 2 z 2 + ... , 
v 

and for density, 

pV N „ N 2 s 

- — = l+B 2 — + 0( — ) + 
NkT V V 



(9) 
(10) 

(11) 



Inspection of Eqs. <j4j and (0 reveals the relation 
B 2 = -b 2 , 



(12) 



While we derived the relation for the ideal Ferm i gas, E q.[T2lis 
a general result, valid for all real gases (see, e.s. iHuana [1963). 
The second virial coefficient is a function of temperature only; an 
explicit expression is obtained from classical statistical mechan- 
ics by the two-body configuration integral over the interaction 
potential U\ 2 



B 2 
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(13) 



Of course, as infinite series, activity and density expansions 
would be strictly equivalent, but when dealing with recombina- 
tion, a density expansion would have to go to a very high order 
to match the accuracy of even a second-order activity expan- 
sion. Indeed, for recombinations, the activity expansion ( 19110b is 
vastly superior to the virial expansion ( fTTT i. This is easily seen 
from Eq. f!3) . where for a (necessarily) attractive two-body po- 
tential and low temperatures B 2 becomes large negative, which 
will cause inevitably a negative total pressure at some low tem- 
perature. In other words, the virial expansion breaks down. This 
breakdown is illustr ated by the dashed line of Figure Q] (from 
Ebeling e t al.[ Il976l) . For simplicity, that figure is for hydro- 
gen molecular dissociation, but qualitatively it would be the 
same for ionization. In contrast, the activity expansion can avoid 
this breakdown, as revealed by a more detailed discussion of 
Eqs. d9U0|l that can be found elsewhere ( e.g. Ebeling et all 
119761: iKremp et al\. 120051) . Importantly, already to lowest non- 
ideal order, the activity expansion can correctly describe com- 
plete low-temperature recombination (solid line in Figure [TJ. 
This property of activity expansions is the principal reason for 
the remarkable success of the OPAL equation of state. 

In contrast, as a virial expansion, the FK equation of state 
suffers from the pathology illustrated by Figure [1] irrespec- 
tive of the fact that its low-order virial coefficients are exact. 
Nevertheless there is a domain of applicability where we can use 
FK: Figure Q] shows the overlap region of the virial and activity 
expansions. A natural quantitative limit for the application of a 
virial expansion becomes apparent: as long as the physical con- 
ditions are such that more than about 90% of the bound particles 
remain ionized, the resulting thermodynamic quantities still will 
be reliable. In the following we adopt this criterion for the FK 
approach. 

3. The Path-integral formalism in the framework of 
the Feynman-Kac representation. 

In many situations, stellar matter can be adequately described in 
terms of quantum plasmas made of electrons and nuclei. In this 
part, we summarize the main steps of the path integral formal- 
ism in the FK representation. Details of the calculation can b e 
found in dAlastuev & Perezl Il992t lAlastuev et ali 11994 119951) . 
The general method can be divided in two steps. In a first step 
only Maxwell-Boltzmann (MB) statistic is considered and all 
exchange effects are omitted. Using the FK representation, an 
equivalent classical system made of closed filaments interact- 
ing via two-body forces is introduced. Applying traditional tech- 
niques of classical mechani cs inspired by Abe-Meeron's method 
(Abe, 1959; M eerorlll958l) . a formal diagrammatic representa- 
tion of the MB quantities was obtained, which is term-by-term 
well-behaved. 
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Fig. 1. Illustration of the behavior of activity (solid line) and 
density expansions (dashed line) for reacting systems, here 
for the simple reaction of neutral hydrogen H recombining 
to hydrogen molecules H2. Ad refers to fully dissociated 
hydrogen. The activity expansion (a) correctly describes re- 
combination with just one non-ideal term (equations 191 1 0b . 
Without further high-order terms, the corresponding virial ex- 
pansion (b) (/.<?. the reacting analog to the single-species equation 
ITTb ex hibits catastrophic n on-physical behavior (negative pres- 
sure!) (Ebeling et alW\9T(h . The case of ionization is analogous. 



In a second step of this approach, the exchange contribution 
was introduced, evaluated perturbatively via Slater sums. Use of 
the FK representation leads to the appearance of open filaments 
(impurities) immersed in a bath of closed filaments described by 
MB statistics. This problem can be dealt with using the method 
that is explained in the MB case to inhomogeneous situations. 

Collecting all the contributions arising from both the MB 
and the exchange treatment, a systematic virial expansion of 
the equation of state was obtained, which includes screening, 
diffraction, bound states and exchange contributions to the ideal 
gas (depending on the statistics of the particles). The expres- 
sion of the pressure was calculated up to and including order 
p 5 / 2 (with p being the overall densit y of the system). To thi s 
order of densi t y, the calcu lation of lAlastuev & Perez] (|1992|) ; 
lAlastuev et al.\ d 1994 1 19951) describes exactly all the effects per- 
taining to bound states, diffusive and exchange processes of a 
system consisting of two species of elementary particles. 
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e a p a ) , A a p = (pti z /m a/j ) , m a p = 

a 

m a mp/(m a + nip), x a p = -y/2l a p/A a p, l a p = /3e a ep, C\ = 
15.205 ± .001, C 2 = -14.733 ± .001, and the Euler-Mascheroni 
constant C = 0.577216. 

In Equation (fT4l . Qix a p) refers to the so-called quantum 
second- vi rial coefficient first introduced by Eb eling and co- 
workers (Ebeling eTall 1 1 976t iKraeft et all [l986h . 
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while E(x aa ) is defined as the two-body exchange integral, 
E(x aa ) = (2 Vtt) lim \ dr < -r\e~ phm \r > (16) 

R ^°° Jr<R 



In Equations ( fTBl and ( fT~6b , h a p is the one-body Hamiltonian 
of the relative particle with mass m a p submitted to the Coulomb 
potential e a ep/ r. The functions Q and E only depend on the tem- 
perature via the single dimensionless parameter (-1/A). 



3.1. Exact expression of the equation of state to the order of 3 2 Derjvjng 0{her Thermodynamic Variables 



The virial expansion o f the pressure to order p 5 ^ 2 gives 
(lAlastuev & Perez! 1 1992I) 
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The expression for the pressure in Equation < TT~4T > can be used 
to derive other thermodynamic variables. For example, the rela- 
tionship between the pressure and the free energy is given by 



op V 3( Pf) Of 

pp = Pf 

V 8 P« 



(17) 



where / is the free energy per unit volume. The internal energy 
per unit volume u is then found as a derivative of the free energy 



u = — — (Bf) a 
dp w p 



(18) 



Other quantities of interest can be derived in a similar way via 
the usual thermodynamic derivatives. 
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Fig. 2. Coupling parameters for electrons, protons, and helium Fig. 5. Degeneracy parameter for helium nuclei from the surface 
nuclei from the surface to the center of the Sun. to the center of the Sun. 
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Fig. 3. Degeneracy parameter for electrons from the surface to Fig. 6. Ionization degree of hydrogen from the surface to the cen- 
the center of the Sun. ter of the Sun. 
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Fig. 4. Degeneracy parameter for protons from the surface to the Fig 7 _ First ionization degree of helium from the surface t0 me 
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3.3. The general structure of /3P 

First, the classical contributions of the long-range part of the 
interactions are polynomials in the inverse temperature, the 
charges and the densities, which do not involve Planck's con- 
stant. The involved coefficients are evaluated analytically or, like 
C\ and C2, by numerical computations of dimensionless inte- 



grals. The lower-order contribution is nothing but the familiar 
Debye-Hiickel term in p 3 ^ 2 , and constitutes the leading correc- 
tion to the MB ideal pressure in Equation ( PT7l >. It can be quoted, 
here, that the classical OCP limit up to order p 5 ^ 2 included does 
coincide with those calculated bv lCohe n & Murphv ( 1969), as it 
should. 
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Fig. 8. Second ionization degree of helium from the surface to 
the center of the Sun. 
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Fig. 9. First and second ionization degree of oxygen (as a rep- 
resentative heavy element) from the surface to the center of the 
Sun. 

The contribution of quantum diffraction at large distances 

appears only at the order p 5 ^ 2 and reduces to Y, a ~ z ~ K £,P a 
The occurrence of this term shows that the long-range part of 
the interactions cannot be entirely treated at a classical level. 
The diffraction correction is merely proportional to Ti, because 
it arises from large distances where the quantum effects can be 
treated perturbatively "a la-Wigner-Kirkwood". 

The term p a pp Q(- V^y 2 -) is the total contribution from 
both bound and scattering states of two charges e a and ep. The 
truncation of < r | exp (-/3h a p) \ r > in the integral defining Q en- 
sures that this contribution is finite. This regularization is not 
an arbitrary mathematical artifact. It is directly related to the 
truncated structure of the bond fj, and reflects the screening of 
the Coulomb interaction at large distances. For opposite charges 
such that e a ep < 0, one may extract from Q a contribution of the 
bound states which reduces to the familiar Planck-Larkin (PL) 
sum 



CO 

2« 2 [exp(-y6e^)-l+ySe„ a/ 



(19) 



tions of the bound states contributions can be introduced from 
Equation ( Tl3T > by us ing the basic properties of the trace. For in- 
stance, as shown by iBollel dl987l) . there exists an infinite set of 
arbitrary decompositions in terms of bound and scattering con- 
tributions of the PL sum itself. So, as far as thermodynamic 
quantities are concerned, only the total contribution of both 
bound and scattering states is an unambiguous quantity. The 
p 5/,2 -contribution from bound and scattering states merely re- 
duces to its p 2 counterpart multiplied by f3e a epKr,. This multi- 
plicative factor arises from many-body effects which induce a 
constant shift -e a epKo on the energy levels of the two-particles 
states. _ 

Finally, the contribution p 2 .E(- V2 j 2 -) arises from the ex- 
change of two charges e a in the vacuum. It is finite, indepen- 
dently of any screening effect, because the off-diagonal matrix 
elements < -r | exp {-f3h aa ) \ r > are short-ranged. The mag- 
nitude of this contribution is smaller than the one relative to 
free particles because the repulsive potential <? 2 /r inhibits the ex- 
change. Similarly to what happens for the contributions of bound 
and scattering states, at the order p 5 ^ 2 , the many-body effects on 
the two-particle exchange amounts to lowering the repulsive bar- 
rier e\jr by the constant -e 2 kq. 



3.4. The nature of the nonideal terms in the FK equation of 
state 



where e"f = -e 2 e^m a p/(2H z n z ) are the energy levels of the 
hydrogenic atom with Hamiltonian h a p. However, other defini- 



As is usual for a weakly nonideal plasma (see, e.g. lEbeling et ali 
1 1976b . the most important contributions beyond the ideal gas 
limit are the exchange terms, the Debye screening corrections, 
the short-range corrections and diffraction (in decreasing order 
of importance). Let us recall briefly the various conditions al- 
lowing us to check the nature of the short-range corrections for 
a given temperature and density. 

Since the current level of our FK equation of state is trun- 
cated at the order p 5 ^ 2 , it can only include processes with two- 
body effects. It comes naturally, from the structure of the equa- 
tion of state, that short-range quantum effects are virtually de- 
coupled from the other corrections. This means that in all orders 
of density, the main impact of these short-range quantum correc- 
tions merely pertains to the recombination and ionization pro- 
cesses. The magnitude of the neglected higher-order terms will 
indicate the influence of three-body or higher interactions. The 
comparison with solar data will allow us a quantitative estimate 
of these neglected terms. 

Similarly, one has to be aware of the higher Debye correc- 
tions beyond the p 5 ^ 2 term. A careful analysis of the virial ex- 
pansion shows that the quantum short-range correction that only 
involves the two-body virial coefficient Q(x a p) results in an al- 
ternating series in p, which is therefore convergent. Since the 
leading Debye correction to pressure (included in our formal- 
ism) is negative, the following term must be positive. Should at 
some physical conditions the neglected higher corrections lead 
to a lowering of pressure, we will be able to interpret them as a 
manifestation of a recombination process occurring under these 
thermodynamical conditions. 

4. Limits of validity of the FK equation of state 

4.1. Dimensionless parameters 

A detailed analysis of this equation of state suggests rewriting 
the truncated expression of the thermodynamic quantities (for 
example pressure) in terms of three characteristic lengths: 
• mean interparticle distance d a = (3/47rn Q ,) 1 ^ 3 
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• thermal de Broglie wavelength A a = (J3h 2 /m a ) 1 ^ 2 

• Landau length l a p - fie a ep 

Here the notation is standard, and there should be no confusion 
between ji being both a label and l/(kT). Note that since we are 
now dealing with quantitative applications, we denote number 
density by the usual symbol n rather than the p that was used in 
the preceding theoretical sections of this paper. 

Before the FK equation of state can be applied to the Sun, we 
must determine where the formalism is applicable. Here we eval- 
uate the relevant criteria throughout the Sun, based on a standard 
so lar model (Model S from Chr i stense n-Dalsgaard, as described 
in lChristensen- Dalsgaard et aTL 1 1 9961) First, we verify the well- 
known result that the plasma in the solar interior is indeed not 
strongly coupled. Figure [2] shows the coupling parameter 

u a 

through the solar model. The index a denotes electrons, protons, 
and helium nuclei, respectively. Nowhere does the coupling pa- 
rameter T a come near to the critical value of unity. 

Figures [3][5] show the degeneracy parameter n a h? a , again for 
electrons, protons, and helium nuclei, respectively. 

4.2. Qualitative discussion of ionization and recombination 

As we discussed in detail in section 2, avoiding too much recom- 
bination is the most restrictive criterion for the FK equation of 
state. Figures |6]|9] show that it is certainly more restrictive than 
coupling and degeneracy. Obviously, our low-density virial ex- 
pansion cannot be used in the cooler outer layers of the Sun, 
where the plasma is far from fully ionized. To obtain an estimate 
of limit of validity of FK, which we identified with 90% ion- 
ization, we have calculated the hydrogen and helium ionization 
fraction using the standard (ground-state-only) Saha equation 



N z+ iN e - 
N z 



gz 



XT' 



throughout the same solar model. As usual, N a and g de- 
note the number of particles and statistical weights of the var- 
ious species (here identified by the ion charge Z) and xz is the 
ionization potential pertaining to the ionization reaction Z — > 
Z + 1. Using this simple Saha equation for this purpose is al- 
right since we use it only for a kind of error analysis. And, as 
numerous studies have shown, for solar conditions, this sim- 
ple Saha equation is everywhere at m ost 10% wrong ( e.g., 
I Christensen-Da lsgaard & Dappenl [l992l) . 

Figures |6]|9] show the recombination fraction for hydrogen, 
helium, and the first two ionization stages of oxygen (as a typical 
representative heavy element). The 90%-ionization level is also 
shown. Clearly, at some point in the colder outside layers, our 
virial expansion must break down. 

However, in the deeper layers of the Sun (and non-compact 
stars in general), we can verify that the weak-coupling and low- 
degeneracy criteria are fully satisfied and the limitation from the 
"recombination criterion" is totally negligible for the tempera- 
tures and densities of those regions. 

4.3. The treatment of heavier elements 

Adding heavier elements (other than hydrogen and helium) re- 
duces the domain of validity of the FK equation of state to even 
higher temperatures. This is obvious because it is harder to ion- 
ize these elements. Fig. [9] shows this quantitatively for oxygen. 



We therefore adopt the simplest solution which is to ignore the 
heavier elements altogether. Indeed, in the solar application of 
Paper II, we will merely add the pressure of the heavier elements 
to the FK expansion[14] Since in the Sun, the number abundance 
of all elements other than hydrogen and helium is only about 
10~ 3 , adding the heavier elements to FK with an equatio n of state 
such as MHD, which is at least as accurate as 1% ( e.g..lD"Ippen. 
2006) will limit the total error of such a simplified procedure to 
an acceptable order of magnitude of 10~ 5 . 

In addition, knowing the domain of validity also allows us 
to study selected recombination reactions in the framework of 
the FK equation of state. For instance, in the calculation of nu- 
clear reaction cross sections, the ionization degree of the various 
nuclei plays a role, because the statistical presence of electrons 
around the nucleus modulates the reaction rate (enhancement 
factor). This has consequences for solar physics especially with 
respect to the neutrino problem, where the r eactions involvin 
Be an d B have to be know very accurately (iPerez & Ch abriei 
2004). 



4.4. Solar models with a hybrid equation of state 

Because of its limitations to nearly full ionization, it is impos- 
sible to use the FK formalism for the entire Sun. Therefore, 
for our applications in Paper II, we will have to use a hybrid 
two-zone equation of state, consisting of FK serving in the inte- 
rior, and a more conventional one (OPAL or MHD) in the outer 
layers. Although such a hybrid equation of state will have dis- 
continuities at the transition, these discontinuities will not pre- 
clude the solution of the usual stellar structure equations. Their 
main effect is expected to be a small jump in the gradient of 
the physical quantities, similar to that caused, for instance, by 
standard interpol ation in equation-of-state and opacity tabl es. As 
detailed studies dChristensen-Dalsgaard & Dapped, [1992) have 
shown, successful equation of state tests can still be made despite 
such small imperfections. The analysis of Paper II will reveal if 
the small discontinuities will cause problems. If so, they can be 
avoided by the introduction of suitable smooth ramp functions at 
the juncture line of the two-zone equation of state. 

4.5. Distribution of FK equation of state 

OPAL and MHD have been made available to solar and stel- 
lar modelers through the use of tables. Although this saves the 
modelers a bit of computing time, the tables have only been pro- 
vided for a limited range of conditions. In order to make the FK 
equation of state more versatile, we plan to develop a computa- 
tional routine for which users can define their own input param- 
eters such as element abundances. This implementation of FK 
will be ideally suited for exploring issues such as He diffusion 
and the inconsistency between seismology and the new values 
of heavy element abundan ces (lAsplundl 120051: lAsplund et ali 
l2005HGrevesse et a/.Ll2007h . 



5. Conclusions 

Despite the limitation to nearly fully-ionized matter, the FK ap- 
proach can be applied to large parts of the solar interior, more 
precisely to the part where helium is at least 90% ionized. In 
those regions the FK approach is applicable, and it is especially 
suited to study screening effects, bound states, recombination of 
fully ionized ions and exchange effects. 
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Paper II will be dedicated to these observational 
tests. Since the important relat ivistic correction for elec- 
trons ( El liott & Kosovichevi 1 19981) is not part of the original 
FK equation of state [14] we will add it in Paper II. Paper III 
will go further, on a rigorous theoretical level, to the calculation 
of the thermodynamic functions up to order p 3 in density. The 
additional higher-order terms will allow us to take into account 
the three-body effects that occur in a moderately ionized plasma 
(for example, in the case of the solar center, the possibility of a 
recombination of He ++ is still an unsolved problem). 

Despite the severe technical difficulties involving the practi- 
cal application of the FK equation of state for solar physics, there 
are unique features in the FK approach that promise to turn it into 
the most exact of the available formalisms, if they are limited to 
the layers of the Sun where more than 90% of helium is ionized. 
The localizing power of helioseismology allows us to test FK 
observationally, promising both better and tighter constraints on 
the equation of state. For instance, by virtue of its exactness, the 
FK equation of state becomes a promising tool for tweaking pa- 
rameterizable equations of state such as MHD. More generally, 
the FK equation of state will lead to better solar models as well 
as tighter solar constraints of the equation of state. 
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